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Abstract 



In this paper, the Slater integrals for a screened Coulomb interaction of the 
the Yukawa form are calculated and by fitting the Thomas-Fermi wavevector, 
good agreement is obtained with experiment for the multiplet spectra of Pr^+ 
and C/^"*" ions. Moreover, a predicted multiplet spectrum for the heavy fermion 
superconductor UPt^ is shown with a calculated Coulomb U of 1.6 eV. These 
effective Coulomb interactions, which are quite simple to calculate, should be 
useful inputs to further many-body calculations in correlated electron metals. 
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Calculations of effective Coulomb interactions are important not only in atomic spec- 
troscopy but also in condensed matter physics where the problem of correlated electron 
metals has received renewed attention with the discovery of heavy fermion and high temper- 
ature superconductivity. Even for an isolated atom, these calculations are difficult given the 
large number of possible configurations which must be considered and in the case of heavy 
ions, the added complication of taking into account relativistic effects. Only recently have 
detailed configuration interaction (CI) calculations been performed on Pr^~^, the proto- 
type ion usually considered when trying to understand f electron correlations. p[ No similar 
calculations exist for f/'^"'", the actinide analogue of Pr^+. The latter ion is of significance 
in condensed matter physics since uranium heavy fermion metals are lattices of such ions. 
1^ In fact, recent studies of the /^-/^ Anderson lattice model show that atomic correlations 
can have a profound impact on the mean-field quasiparticle bands. 0] Moreover, it has been 
proposed that such correlations also appear in the residual quasiparticle-quasiparticle inter- 
actions and are responsible for heavy fermion superconductivity. If anything, a proper 
understanding of these atomic correlations will be necessary to achieve a complete picture of 
heavy fermion physics and probably copper oxides as well. In this paper, a simple method 
for calculating these correlation effects is presented. The screened interaction is taken to 
be of the Yukawa form and is evaluated using relativistic wavefunctions from a Hartree- 
Fock program. The calculations are then fit with a single parameter (the Thomas-Fermi 
wavevector) to the observed multiplet spectra of Pr^+ and U^^ . Moreover, results of recent 
high energy neutron scattering data are used to make a prediction of the multiplet spectrum 
for the heavy fermion superconductor UPt^. 

The standard method for interpreting multiplet spectra in heavy ions is to perform a 
least squares fit of the spectra using a Hamiltonian with adjustable coefficients. The 
main terms for /" spectra are the effective Slater integrals (L=0,2,4,6) and the spin- 
orbit splitting. Smaller terms corresponding to odd L due to configuration interaction plus 
other magnetic terms, which are necessary to get an exact fit to the spectra, are ignored here. 
These effective Slater integrals are considerably smaller than their Hartree-Fock values. This 
is as expected since the effect of correlations is to screen these values (in CI language, the / 
electron is distributed in other configurations, thus reducing the wavefunction overlap and 
decreasing the value of these integrals). It is difficult, though, to calculate these effective 
integrals from first principles. Morrison and Rajnak have shown that the corrections from 
many-body perturbation theory have a slow convergence for Pr^^ . Recently, Cai et al 
have performed detailed configuration interaction calculations for Pr^+ including up to 1708 
configurations for each J state. The agreement with experiment is significantly improved 
over Hartree-Fock, but there are still some discrepancies. 

Here, an approach intermediate between the fitting procedure and that of the ab initio 
calculations is taken. The reason for the problems mentioned above is that both ab initio 
methods are rather inefficient in calculating these effects. In solid state physics, one has to 
sum infinite series to describe screening effects properly, so the lack of convergence at third 
order in the above many-body perturbation calculation is not too surprising for an atom 
with a large number of electrons like Pr. As for CI, wavefunction expansions of this sort 
have very slow convergence except for small atoms when most of the correlation is coming 
from near degeneracy effects such as Be. Again, these problems are more severe when large 
numbers of electrons are involved. The most efficient way of dealing with this, then, is to 
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assume a screened interaction from the outset. In particular, in this paper, a Yukawa form 
is used 

g-A|ri-r2l 

V{n,f2) = — — — (1) 

|ri - r2| 

The work of this paper is still semi-phenomenological in that A, the Thomas-Fermi wavevec- 
tor, is treated as an adjustable constant. Ab initio calculations of it would suffer the same 
slow convergence problems as mentioned above for the many-body perturbation calculation, 
but at least our ignorance has been transfered to just one parameter as opposed to sev- 
eral of them as in the least squares fitting to the experimental spectra discussed above. In 
solids, reducing the fitting parameter to one number is highly desirable since the detailed 
information available in atomic spectroscopy is often lost due to hybridization effects. 

The Slater integrals of the Yukawa model are easily found. Just as the bare Coulomb 
interaction is expandable in powers of r, the screened interaction is expandable in spherical 
functions 

y(fi,f2) = -A^(2L + l)ji(2Ar<)/i?)(ar>)Pi(cos^) (2) 

L 

where ji is a spherical Bessel function, h^l^ is a spherical Hankel function of the first kind, 
and Pl is a Legendre polynomial, with 9 the angle between fi and r2 and r< (r>) the 
lesser (greater) of the two. These spherical functions were determined using routines in 
Numerical Recipes. The wavefunctions used in the Slater integrals were determined 
from relativistic Hartree-Fock calculations. For convenience, a single integral is defined 
by taking a weighted contribution of the two spin-orbit wavefunctions, that is 0^ in the 
non-relativistic Slater integrals is replaced by |05/2 + f07/2- 

In Fig. 1, the evolution of these Slater integrals versus A for the J/^"*" case is shown (A=0 
is the Hartree-Fock result). The L=0 integral is strongly screened whereas the higher L 
integrals are less affected. This is to be expected since the former is a charge fluctuation 
integral whereas the latter are shape fluctuation integrals. For all L, the larger the L, the 
smaller the screening effect, as expected. 

The matrix elements of the secular matrix for the p ion case are listed by Goldschmidt. 



Ill There are 13 eigenvalues (3 triplets, 4 singlets, each of the triplets being split by 
spin-orbit). The spin-orbit splitting determined from the relativistic Hartree-Fock orbitals 
matches the experimental estimate for Pr^+ to within 0.1%. This is as expected since 
the spin-orbit interaction is a one-body force. Implications from the fits that this interaction 
is screened was incorrect because non-relativistic estimates had been used for the Hartree- 
Fock value. A slight discrepancy was found in the [/'^+ case in that the Hartree-Fock value 
underestimated the experimental fit by about 2%. This is probably due to the strong in- 
termediate coupling (LS versus jj) nature of the U ion. For the purposes of this paper, the 
calculated Hartree-Fock values of ^ are used (where 7/2 ^ is equal to the splitting of the 5/2 
and 7/2 levels) as opposed to the "experimental" ones. These values are 0.0990 eV for Pr^^ 
and 0.2387 eV for U^^ . The secular matrix is then diagonalized using Slater integrals for 
a particular value of A. The energy of all levels will be refered to that of the ground state 
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Results are shown for Pr^"*" in Fig. 2 with Slater integrals listed in Table 1. An RMS error 
of 82 meV is obtained for A=2.0. If the highest level (^^So) is ignored, an RMS error of 40 
meV is obtained for A=1.9. This is to be compared to the RMS error of the CI calculation 
|Jl[| which is 163 meV. An interesting point is that the ordering of the levels agrees with 
experiment, therefore questions raised about the assignment of the ^Iq level are probably 
not valid. Another point is that the values of the L=2 and L=4 Slater integrals agree with 
the fitted (experimental) values. On the other hand, the L=6 one is substantially larger. 
The quoted experimental value may be somewhat misleading, as it is more screened over 
Hartree-Fock than than the L=4 one, whereas one would conceptually expect that higher L 
integrals are less screened than lower L ones (as is indeed reproduced here). This leads to 
the possibility that there is another least squares fit to the data which more closely matches 
the current predictions. Such a fit would presumably have different values of the odd L 
interactions than listed in the published fits. This issue will be discussed in more detail 
in the next paragraph. Finally, the calculated L=0 Slater integral (the Coulomb repulsion. 



U) is close to the experimental estimates of 5.3 eV [T^ and 5.5 eV [0. This is quite 
interesting since the fit of A does not involve U. This suggests that the current method 
provides a trivial and independent way of estimating the Coulomb U. These estimates of U 
are important inputs to many-body calculations in solids. 

In Fig. 3 and Table 1, results are shown for f/'^"'". For A=1.6, an RMS error of 84 meV 
is obtained. If the highest level (^5*0) is ignored, an RMS error of 40 meV is obtained for 
A=1.5. These deviations are identical to that found for Pr^+, with A smaller in this case 
since the 5f orbitals are more delocalized than the 4f ones. The level ordering agrees with 
experiment except for an interchange of and ^G^. The level ordering in this energy range 
(including ^Pq) is rather sensitive to A, and this minor problem is almost certainly due to 
the intermediate coupling nature of this ion. This has been checked by directly replacing the 
Slater integral routines in the relativistic Hartree-Fock code by the Yukawa ones (there is a 
large number of Slater integrals due to the various combinations of the two spin-orbit radial 
functions). In this case, the ordering of the two levels comes out correct, but at the expense 
of having the ^Pq lower than these two, as opposed to higher as in experiment. In fact, the 
RMS errors are larger with this approach; in particular, the splittings within an LS term are 
worse than with the other approach (these splittings have a weak dependence on A). This 
again points to the intermediate coupling nature of the ion as the problem. Finally, as in 
Pr^~^, the predicted L=6 Slater integral is much less strongly screened than implied by the 
experimental least squares fits. |]T^ On the other hand, the value reported by Goldschmidt 
|TT| is 0.3 eV larger than that quoted by Van Deurzen et al. This difference is due to 



different definitions of the odd L interactions which in turn alters the values of the even L 
Slater integrals. M Inclusion of these odd L terms improves the RMS error of the current 



work by over a factor of 3, not surprising since 3 more fit parameters have been added. ||T5 
First principles calculations of these odd L terms have had mixed success, [0 so their use is 
somewhat dangerous. 

Also reported in Fig. 3 and Table 1 are results for A=2.6. This value was chosen to fit 
the ^H^ - splitting in UPt^ observed with high energy neutron scattering. [T^ This can 
then be taken as a prediction for the entire multiplet spectrum in this case. It should be 
remarked, though, that given the itinerancy of the f electrons in UPt^, such multiplet effects 
will be very weak. They will probably be better observable in UPd^, where the f electrons 
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are known to be localized and which has an almost identical observed splitting. Of more 
interest is the U value, which is necessary in many-body calculations. The predicted free 
ion value of 3.3 eV is larger than the experimental values quoted by Brewer |12[ of 2.3 and 
2.6 eV, although it should be remarked that these values were extracted from neutral ions 
with an /'^ ground state, and thus this U is expected to be smaller than the one appropriate 
for the more localized U^^ ion. Also, U is predicted here to be reduced by a factor of two 
in the metallic environment appropriate for UPt^ and UPd^. This U value of 1.6 eV should 



be a useful input for many-body calculations in uranium heavy fermion metals. [|T7 

In conclusion, a simple method has been proposed for determining the effective Coulomb 
integrals needed for many-body calculatons, including the Coulomb U. This method has 
been applied to Pr^+ and t/^^ ions and gives reasonable values for these parameters. In 
turn, this method was used to give a prediction of such values for the case of uranium heavy 
fermion metals. 
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TABLES 



TABLE L Slater integrals (eV) for Pr3+, C/'^+, and UPts. In column 2, HF is Hartree-Fock, 
exp is experiment, and the numbers refer to values of A. The experimental estimates were taken 
from Ref. [|l2|,|| {Pr, L=0), Ref. @ {Pr^+, L>0), Ref. {U, L=0), and Ref. [|l| {U^+ , L>0). 







po 


F2 


i74 


F6 




HF 


25.722 


12.227 


7.670 


5.517 




L9 


5.560 


9.000 


6.811 


5.190 




2.0 


5.230 


8.784 


6.735 


5.158 




exp 


5.3,5.5 


9.090 


6.927 


4.756 




HF 


18.796 


9.482 


6.202 


4.552 




1.6 


3.309 


6.377 


5.281 


4.185 




exp 


2.3,2.6 


6.440 


5.296 


3.441 


UPh 


2.6 


1.645 


4.421 


4.336 


3.724 
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FIGURES 



FIG. 1. Slater integrals (eV) versus A (Thomas-Fermi wavevector in a.u.) for the U'^'^ ion. 

FIG. 2. Multiplet spectra (eV) relative to the ^H^ ground state of Pr^~^ for Hartree-Fock (HF), 
A=2.0 (2.0), experiment (exp) and configuration-interaction (CI) [Q]. The first excited state is 

FIG. 3. Multiplet spectra (eV) relative to the ^H^ ground state of C/^^ for Hartree-Fock (HF), 
A=1.6 (1.6), experiment (exp) [O], and A=2.6 (2.6). The latter A was fit to the observed multiplet 



splitting in UPt^. |16] The first excited state is F2 except for Hartee-Fock where it is H^. 
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